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Quark contact interactions are an important signal of new physics. We introduce a model in which 
the presence of a symmetry protects these new interactions from giving large corrections in flavor 
changing processes at low energies. This minimal model provides the basic set of operators which 
must be considered to contribute to the high-energy processes. To discuss their experimental signa- 
ture in jet pairs produced in proton-proton coUlisions, we simplify the number of possible operators 
down to two. We show (for a representative integrated luminosity of 200 pb~^ at ^/s = 7 TeV) how 
the presence of two operators significantly modifies the bound on the characteristic energy scale of 
the contact interactions which is obtained by keeping a single operator. 

PACS numbers: 13.85.Hd,12.60.Rc., 14.65. Jk 



I. MOTIVATIONS 

Fermions like quarks can be made to interact directly — that is, without the exchange of an intermediate particle — by 
simply adding to the standard model (SM) lagrangian four-quark contact terms like, for example, 

27T - 

IpLl'^lpL'il^Lli.i^L , (1) 

where 27r gives the strength, A the characteristic energy scale of this new interaction and the quark fields ip are taken 
to be left-handed. 

Because of the non-renormalizability of the term in eq. ([T]), such an operator is often thought as the low-energy 
effective approximation of a renormalizable lagrangian in which heavy particles are exchanged. These heavy particles 
can be of many different kinds, each kind giving rise to an effective operator with different color, flavor and Dirac 
structure. Sometimes these heavy states are thought of as a substructure of the quarks themselves and in this case the 
contact interactions are presented as evidence for quark compositeness. More in general, the heavy states represent 
new physics which lives at an energy scale that is too high to manifest itself with the production of the new states 
either as intermediate resonances or in chain-decay processes, and the effect of which can only be seen by the effective 
operators of the contact interactions [Tj . Non-renormalizability is not a problem in models in which the couplings run 
toward a ultraviolet fixed point. In these asymptotically safe (AS) models [2] the contact interactions in eq. ([T]) can 
be considered as fundamental. It has recently been shown that indeed such operators arise in a natural manner in AS 
models of the weak interactions and a search for their presence could provide an important experimental clue [3]. 

In any case, be the contact interactions fundamental or remnant of new physics at higher energies, the search for 
their existence and the bounds on their characteristic energy scale is important. The LHC has already provided us 
with new constraints [11 [5] . These constraints were derived by assuming the existence of only one kind of contact 
term, namely that in eq. ([ij in which the operator is given by the product of two left-handed quark currents. This has 
become the standard practice following [I] because it is simple. Unfortunately, this is too restrictive an assumption 
and the significance of an analysis based on it is unavoidably weakened. 

It is easy to imagine a great variety of different operators contributing to the quark contact interactions. The 
problem is that this variety is constrained by stringent bounds on flavor physics at low energies [B]. We do not want 
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to track down every and each operator for its possible low-energy effect and, in order to provide a minimal model, 
we impose a symmetry on the possible contact interactions which makes them safe with respect to these low-energy 
constraints. This model defines a basic set of operators whose size is not severely constrained by flavor physics and 
the existence of which can be tested in high-energy processes. It is also general enough to make the bound on the 
characteristic energy scale A realistic. 

In principle, all these operators should be entered in the analysis, each with a different strength. To make the 
analysis manageable, we simplify further down to two the number of possible operators. This skeleton model is 
sufficient in showing — in jet pairs produced in proton-proton colllisions at the LHC for a representative integrated 
luminosity of 1 fb~^ at i/s = 7 TeV — how the bound on the energy scale of the contact interactions is sensitive 
to the relative strengths of different terms. This exercise shows that at least a subset of operators should be taken 
into account, current analyses based on a single operator cannot be considered as final, and the actual bound on the 
characteristic energy scale A is weaker than reported [U [5] . 



II. THE MINIMAL MODEL 

Let us for a moment consider the case of a single fermion family with quarks with the same mass. The most general 
four-quark interaction will depend on the overall symmetry we want to impose on the system. It could be U{1)q or 
SU{2) X U{1) and in these cases we would have, respectively, 20 or 10 possible terms. Following the general idea that 
stronger interactions are more symmetric than weaker ones, we want to be more restrictive and impose a larger group 
SU{2)l X S'[/(2)fl and a parity symmetry. 

The complete set of SU{2)l x SU{2)ji and parity invariant four-fermion operators is given by four independent 
terms: 

+ A4 (^r7M^Lv;iVvr+vik°7MV'k'vi^jV^^r) (2) 

where the dimensional coefhcients A„ can be written as 27r/A^j in terms of the characteristic energies A„. Indices i,j 
and a, h are SU{2) and color indices, respectively. 

The operators in eq. ^ must be considered together with the Yukawa term of the standard model: 

-j{i^^U^'^^^+h.c). (3) 

For a realistic model wc should consider the splitting between up and down-type quarks as well as the three SM 
families. This gives rise to a large proliferation of possible terms. At the same time, we must take care that the 
four-fermion operators do not yield unwanted flavor-changing neutral current (FCNC) processes with /S.F = 2 which 
are strongly suppressed by the experimental data. 

In such a realistic and most general case the A^ in eq. ^ are 4 flavor-index tensors and the coefficient h in eq. ([s]) 
is a matrix. Thus the four-fermion interaction lagrangian has 20 4-index tensor operators, that in general are not 
simultaneously diagonalized with the Yukawa mass terms, thus giving rise to problematic operators such as 

{dLSRf , (4) 

which affect for example meson oscillations [5]. 

To prevent such operators and reduce the number of free parameters one possibility is introducing a flavour sym- 
metry. In the following we assume that left and right handed quarks transform as the fundamental representation of 
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a continuous family symmetry we choose to be SU{3)f- In this case the set of four fermion operators given in eq. ^ 
becomes 



r' 
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jb lia 



(5) 



where a,/3 are SU{3)f flavor indices. The symmetry of the four-fermion interactions is thus SU{2)l x SU{2)ji x 
SU{Z)f. 

In this minimal model there are 8 arbitrary coefficients — the A„ and A„ in eq. ([s]) — each of them multiplying various 
operators which are different for flavor and color structure. In a numerical study, all coefficients should in principle 
be varied and the most relevant among the operators included. 

For what concerns the Yukawa term given in eq. ([s]), first of all we have to split up and down- type quarks: 



2/1" _ 2h'^ 



f 



f 



h.c), 
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where the projectors Pu.d project on the up and down components of ip-'j^ respectively. 

In the most general case eq. ^ breaks the full group SU{2)l x SU{2)i^ x SU{3)f to the electric charge U{1)q in a 
complete arbitrary way. However, if we think the h'^'^ as arising from the vacuum expectation values (VEV) of a field 
Y, they can be written as /i"''^ = (Y'^''^') /Af and, in this case, the Yukawa lagrangian that leads to eq. ^ presents the 
extra accidental global symmetry SU{3)'^^ x SU{3)p^^ x SU{3)p^ according to which Y'^-'^ transform as the (3,3, 1) 
and (3,1,3) representations, respectively. We may now assume that F"''* develops VEV only along the diagonal 
direction SU{3)p^^ x SU{3)p^: in this way, the Yukawa mass matrices h"^ {q — u,d) are symmetric and ~ V^. 
Notice that V^'"^ satisfy Vckm — Vl^Vl^ ^i^*^' thanks to our assumptions, Vckm — V^^V^. Thus in the Yukawa 
sector the full symmetry 5'C/(3)^^ x 5;7(3)^^ x SU{5)j,^ D SU(?,Yp^ x 5'C/(3)f,^, is broken to [U{lYfp x [U{lY]%: 

It is now simple to identify whether a term of the four-fermion operators of eq. ([5| may give rise to flavor violation 
or not: those terms that are invariant under SU{3)p^^ x SU{3)p^ do not violate flavor because are simultaneously 
diagonalized with the Yukawa couplings, the others do. 

Indeed, if we classify quarks according to their electric charge q and their flavor charge qF we see that when in the 
four-fermion operators flavor indices are contracted between quarks of the same electric charge q automatically the 
total qF is zero for each flavor involved. This happens for exactly 14 operators (each with 2 indices running over the 
3 families): 

• 2 operators from Ai(A2) and 4 operators from Ai(A2), 2 of which have the same structure of those from Ai(A2), 
for a total of 8 operators; 

• 2 operators from A3(A4) and 3 operators from A3(A4), 2 of which have the same structure of those from A3(A4), 
for a total of 6 operators. 

These operators correspond to a subset of possible contractions and give rise to the lagrangian 



4"^ = (Ai+Ai)[<,u^> 



La "-Ra '^Re 
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+ (A3 + A3) [ul^l.ul^ul^rul, + dl^l.dljl^^^dl^] 

+ A2 [sL4«4.<. +4.4.4.^2 J +A4s2.7M"Lrf"L7''d2. 

+ (A4 + A4) [ul^lAAul'^^s + 4.7^4.4.7^^^] ' (7) 

where A, B are mass eigenstates indices, while a, 6 color indices (notice the different color contractions in the various 
similar operators). This lagrangian has, as before in eq. ([5| 8 coefficients and f4 operators. Each of these operators 
generates 6 terms once the 3 families are included. 

On the contrary when flavor indices are contracted between quarks that have |Ag| ~ 1 the operators present a total 
(7f 7^ 0. This is the case for 6 operators 

• 2 operators from Ai(A2), for a total of 4 operators; 

• 1 operator from A3(A4) for a total of 2 operators. 



The corresponding operators become 4 flavor-index tensors, the structure dictated by the Vckm entries according to 
\Vckm)ab{vIk,,)cd {ulMi^Rc<n) + {vIkm)ab(Vckm)cd {dlXRs)i^RcdRo) 



- Ai 



+ A3 {Vckm)ad(V^km)bc «,7m<J(4c7^4„) 



+ A2 



{VcKM)ABiVij,M)cD «,4j(4c"flJ + iViKM)AB{VcKM)cD {diyj,J{u'h,^d'i,J 



+ A4(ycKM)AD(V2^M)BC«,7A.<J(4c7^d2J- (8) 

For the operators in eq. ([8| flavor breaking is manifest since they are characterized by qx — ±1 with X ^ A,B, C, D. 
However these operators may mediate FCNC AF = 2 processes only at one loop. Thanks to the suppression related 
to the CKM entries their effect is of the same order, or even smaller, than the SM contributions. 

To summarize: by imposing an appropriate flavor symmetry, the large number of four-fermion operators present in 
the most general case have been grouped in two classes: one that includes operators that conserve flavor, the other 
those that violate it. These operators multiply different combination of the 8 independent parameters A„ and A„ 
in eq. ([s]). Flavor violating operators are modulated by the square of CKM entries and are therefore suppressed. 
When considering pp collisions, the flavor-conserving operators in eq. ([t]) dominate on those in eq. (|8| — which as 
a consequence can be neglected. Moreover, among the operators in eq. ([t]), the largest cross sections are given by 
those involving only the first family in the initial state and the first and, possibly the second family, in the final 
state. Therefore, a complete analysis of this minimal model should take into account at least the 14 operators in the 
lagrangian £|S of eq. ( 7 ) involving the first family to set bounds on the 8 parameters A„ and A„ . 



III. EVENT GENERATION AND ANALYSIS OF THE SKELETON MODEL 



Dijet production in proton-proton collisions pp — > jj+X is the best channel to search for quark contact interactions. 
In QCD, the jet production rate peaks at large rapidity y, because the scattering is dominated by i-channel processes. 
The rapidity is defined as y = ^ ln(£' + pz)/ {E ~ pz), where E is the energy and p^ the z-component of momentum of 
a given particle. On the other hand, quark contact interactions produce a more isotropic angular distribution leading 
to enhanced jet production at smaller values of |?;|. For this reason, searches for contact interactions at the LHC use 
quantities computed from these dijet rapidity distributions in the high invariant dijet mass (mjj) region. 

In current analyses [HIS], the contact interactions are parametrized by a single standard operator, the one introduced 
in [T] — where however a larger set of operators was introduced and the single-operator scenario was only advocated as 
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a simplification. As we argued in the previous section, a standard set of operators can be identified. The simulation of 
the 14 operators and 8 parameters is rather CPU time consuming and here we only consider the effect of the presence 
of more than one operator on dijet production by introducing a skeleton model that admits just two four-fermion 
operators, which we chose to be the first and the third of eq. ([2|. The size of their couplings are parametrized by 
the characteristic energy scales Ai and A3. While a full analysis is certainly necessary, such a radical simplification is 
already sufficient in showing how the presence of more than one operator gives rise to substantial interference effects 
which modify the bounds on the characteristic energy scale. 

The lagrangian of this skeleton model, written in terms of these parameters, reads: 

2-K - - - - 27r - - 

= Cqcd + ^ {i'LJtj.^L^Ll^IpL + i'RHtii'R'^^Rn^'^R) + ^ {lp L^^ R^J l) ■ (9) 

where the isospin and the color contractions, that can be read out directly from eq. (|2|, are omitted. 

This choice of considering just two operators is useful because it reduces the CPU time for the simulation, simplifies 
the analysis, and shows how it is essential to consider more than one operator. 

We use MADGRAPH V.4.5.0 to simulate LHC dijet production in pp collisions at ^/s = 7 TeV, Monte Carlo samples 
are generated for pure QCD and for QCD modified by the new four fermion interaction terms in eq. Since 
MADGRAPH is a leading-order generator and does not support non-renormalizable interactions, we have implemented 
them effectively, introducing a set of fictitious gauge interactions acting only on the first quark family. In this case, 
the identification is ■0 = {u d), and the mass of the fictitious gauge boson has been choosen to be very high 100 
TeV). The generated events are then passed through PCS, the detector simulator in which the parameters are set to 
reproduce the ATLAS detector performance. 

In order to restrict the simulation in the kinematical region of interest, we have applied the following cuts at the 
generator level: 

• mjj > 1000 GeV ; 

• > 30 GeV; 

• |r/J^*| < 2.8 . 

The pseudorapidity is defined as t] = — ln(tan6'/2), where 9 is the angle between the jet and the beam direction in 
the laboratory frame. 

We have generated Monte Carlo samples for different values of the energy scales Ai and A3 between 1 and 10 TeV. 

Given the cuts described above, the Monte Carlo leading order cross section for each choice of the parameters turns 
out to be (7 ~ 5 X 10"^ pb. An integrated luminosity of 1 fb^^ has been generated for each of the points in the (Ai, 
A3) plane described above. 

The variable x is the quantity used for the angular distribution study. It is defined as a function of the rapidities 
of the two highest px jets in the event, yi and y2'- 

X = exp(2|2/*|), (10) 

In the massless particle limit the center-of-mass (CM) rapidity y* — (j/i ~y2)/2 is used to determine the partonic CM 
angle 9* given the relation y* ~ ^In (^x^jcosl*!) ' ^'^'^^her variable used in the following, derived from the rapidities of 
the two jets, is = (j/i +?/2)/2. The x variable is useful for the comparison of angular distribution predicted by new 
processes with those of QCD. In QCD, gluon exchange diagrams have approximately the same angular dependence as 
Rutherford scattering and so the dN/dx distribution is constant in x- On the other hand, the angular distributions 
of new processes are more isotropic, leading to additional events in the low x region. Subleading diagrams in QCD 
can also rise slightly the distribution at low x- 
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The second important kinematic variable is the dijet invariant mass rUjj, which is also the CM energy of the partonic 
system. It is computed from the two jet four-vectors as 



Too 



(11) 



where E and p are the energy and momentum of the jets. 

Events with at least two jets are retained if the highest pt jet satisfies p!^ > 60 GeV and the second highest one 
satisfies p^r^ > 30 GeV. This asymmetric thresholds avoid suppression of events where a third jet has been radiated, 
while 30 GeV threshold ensures that reconstruction is fully efficient for both leading jets. Events with an additional 
poorly measured jet with px > 15 GeV are vetoed to avoid possible incorrect identification of the two leading jets. In 
addition, x distribution are accumulated only for events for which < 1.10 and \y*\ < 1.70. 
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FIG. 1: X distribution, as defined in eq. (10 1 for events with rrijj > 1200 GeV. Empty points represent the x distribution 
for pseudo-data, the solid line for pure QGD while the dotted line for one specific point (Ai = 1 TeV, A3 = 12 TeV) in the 
generated grid. The various distributions are normalized to an integrated luminosity of 200 pb^^. 



The measure of the isotropy in the dijet distribution, introduced in [5], is given by the variable F^. It measures the 
fraction of dijcts produced centrally versus the total number of observed dijets in a specified dijet mass range: 



A^cvcntsdZ/l < 0-6) 

A^cvc„ts(|2/1 < 1-7) ' 



X AT \ ^ -I '7\ ' ^ ' 



where iVcvonts is the number of candidate events within the y* interval. The central region which is expected to be 
most sensitive to new physics is defined by the interval \y*\ < 0.6 and corresponds to x < 3.32, while |y*| < 1.7 
extends the angular range to x < 30. 

The presence of possible contact terms is tested for each value of Ai and A3 in the highest dijet mass bin: rrijj > 1200 
GeV. For a given pair of values of the energy scales, the corresponding value of is obtained starting from the x 
distribution as in Fig. [T] We have generated a QCD Monte Carlo sample, corresponding to an integrated luminosity 
of 200 pb^^, to be used as a pseudo-data sample. In Fig. [l] the x distribution is shown for dijet events passing the 
selection described above with the additional constraint that the invariant mass of the two hardest jets is larger than 
1200 GeV. The pseudo-data sample (empty points) is well described by the QCD (sohd line), while the contribution 
of a contact interaction term (dotted line), corresponding to the point (Ai — 1 TcV, A3 — 12 TeV) in the grid, 
clearly shows a peak in the low x region, giving then a larger values for (defined in eq. (12)) with respect to both 
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FIG. 2; Isotropy of the dijet distribution as given by the variable F-^. The area outside the curve represents values above the 
value computed from the QCD pseudo-data. The area inside represents — at the 95% confidence level — values compatible with 
the pseudo-data obtained for different energy scales Ai and A3 of the contact interactions. The lower bounds for Ai and A3 
are, respectively, 2.2 and 5.1 TeV. 



pseudo-data and QCD samples. A full set of pseudo-experiments, has then been made for each of the points in the 
grid in order to construct one-sided 95% confidence level. 

The result of this analysis is shown by the contour plot of Fig. [2j The value of extracted from the pseudo-data 
represents in Fig. [2] the level below which the contact interactions are compatible with pseudo-data. The values of 
Ai and A3 which satisfy — at the 95% confidence level — this bound are represented by the area inside the curve. By 
inspection, we find that the lower bounds on the contact interaction scale are given by the values Ai ~ 2.2 TeV and 
A3 — 5.1 TeV, respectively. 

The standard one-operator analysis — which corresponds to taking Ai very large — would give a limit A3 = A « 5.6 
TeV which corresponds to the upper margin of the not excluded area in the contour plot in Fig. [2] The bounds we 
find in the case of two operators are weaker than this one because of interference effects. 

If one insists in having a unique energy scale A even in the presence of more operators, it possible to provide it by 
combining the characteristic scales A„ by means of, for instance, the definition 
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(13) 



This definition has the advantage of providing a bound close to the lowest one and to go into the single-operator limit 



when all scales but one are taken to be large. In our skeleton model with just two operators the definition in eq. ( 13 ) 
gives a bound A 2.0 TeV — a value again weaker than what found in the single-operator analysis. 
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